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Abstract 

Taking into account both the Jahn- Teller lattice distortion and the on-site 
electronic potential fluctuations in the orbital-degenerated double-exchange 
model, in which both the core-spin and the lattice distortion are treated clas- 
sically, we investigate theoretically the metal-insulator transition (MIT) in 
manganites by considering the electronic localization effect. An inverse ma- 
trix method is developed for calculation in which we use the inverse of the 
transfer matrix to obtain the localization length. We find that within reason- 
able range of parameters, both the lattice effect and the potential fluctuation 
are responsible to the occurrence of the MIT. The role of the orbital configu- 
ration is also discussed. 

PACS numbers: 71.30.+h, 71.38.-k, 72.15. Gd 
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I. INTRODUCTION 



Colossal magnetoresistance (CMR), an effect discovered in mixed- valence manganites 
Ri-i^Aj^MnOs (where R and A are rare and alkaline-earth ions, respectively) have attracted 
much attention [1-5]. In the hole doping range 0.17 < x < 0.4, such an effect occurs in the 
vicinity of the metal-insulator transition (MIT) temperature Tmi at which the resistivity 
of the system varies drastically. It is generally accepted that the anomalous transport 
phenomena in these manganese systems are closely related to their magnetic properties, in 
particular the paramagnetic (PM) to ferromagnetic (FM) phase transition upon cooling, as 
in most experimental measurements, Tmi and the magnetic Curie temperature Tc are very 
close to each other. 

The double exchange (DE) mechanism [6] is believed to play a major role in the magnetic 
transition, in which the carriers are ferromagnetically coupled to the Mn core spin due to a 
strong Hund's coupling, resulting in the hopping amplitude depending on the relative Mn 
spin orientation. On the other hand, the mechanism responsible for the MIT has not yet been 
very clear [7]. It was found that in single orbital calculations the random hopping from the 
DE mechanism alone is not enough to drive the MIT [8-10]. Many further theoretical efforts 
have been made to understand the MIT. By taking into account the potential fluctuations 
experienced by carriers due to the random distribution of R^^ and A^'^ , Sheng et al. [8] 
proposed that the electrons are localized above Tc due to strong on-site and spin disorder, 
while at low temperature the alignment of the spins reduce the spin disorder and the electrons 
are delocalized. In their calculation, however, the strength of the on-site disorder required 
for the localization, W > 12t (W is the range of the random on-site energy distribution, 
and t is the one-orbital hopping amplitude), is apparently overestimated, as the reasonable 
magnitude of t is roughly between 0.2 and 0.5 eV [11], while the potential fluctuations 
would amount to 1.7 eV if unscreened, and be further reduced when screening is accounted 
[4]. MiUis et al. [12] used the dynamical mean- field method to study the coupling of the 
carriers to local Jahn-Teller ( JT) distortions and to the Mn core spins, but an MIT was found 
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only at half filling {x — 0). In recent years, there are also much interest on phase separation 
scenario [13-15], in which the metallic and insulating regions coexist in the system and the 
transport is a percolation problem, near the percolation critical point, a small change of the 
fraction of the metaUic regions can induce the MIT. But why this fraction change occurs 
with the variation of temperature still needs an explanation. Very recently. Verges et al. [16] 
proposed a lattice-spin mechanism to explain the MIT, by using Monto Carlo simulation in a 
6^ lattice for a single-orbital DE model coupled with phonons, they found that at x — 0.08, 
in a narrow region of the coupling parameter, an MIT coincidence with the FM to PM 
transition upon heating. 

Based on previous experimental and theoretical studies, it appears that the following 
scattering mechanisms are important to the transport in manganites: (i) the strong DE 
interaction between carriers and the localized spin, which is usually regarded as the basic 
mechanism in manganites, (ii) the potential fluctuations experienced by the electrons due 
to the random and ^4^+ ion distribution, according to Ref.[4], this fluctuations might 
be comparable to the electron band-width, and (iii) the JT distortion which was revealed 
in many experimental measurements. Besides the basic DE interaction (i), in Refs.[8] and 
[16], (ii) and (iii) were considered, respectively. Yet both theories have difficulties in their 
explanation of MIT (overestimated W in the former and not at the right doping in the latter) . 
Furthermore, in recent years it was also realized that in manganites the twofold Cg orbital 
degeneracy and the unique Slater-Koster form of electronic hopping [17] are important to 
many aspects of the system [18]. In many theoretical treatments, such as in Refs[8] and 
[16], to simplify the calculation a single orbital model was adopted, as it was believed that 
the JT distortion would split the Cg degeneracy and reduce the problem to an effective 
single orbital one. However, whether the energy splitting A is large enough to ignore the 
degeneracy is questionable, in Ref.[16] the electron- lattice coupling responsible to the desired 
MIT is around A = 1.45, for such a strength of JT couphng, the average of the corresponding 
A is about 4(l-x)t, considerably smaller than the width of the energy band 12t. On the 
other hand, even in the A ^ oo limit where at each site there is only one orbital valid, it 
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is still different from the usual single orbital model due to the anistropy of the hopping, for 
example, in this large A limit if the lower-energy orbit is (1^2 _y2, then along the z-direction 
there is no electronic hopping, very different from the isotropic hopping in the usual single 
orbital model. 

Because of the problems in previous MIT theories raised above, it is desirable to inves- 
tigate the MIT of manganites by including all the three scattering mechanisms as well as 
the orbital degeneracy. Up to now no such theoretical study has been made yet, and we 
will make such an effort in this work. Starting from an effective degenerated DE model 
with electron-lattice interaction and potential fluctuation included, we investigate the MIT 
by considering the electronic localization effect. An inverse matrix method suitable for the 
present specific Hamiltonian is developed for our calculation. We find that while similar 
to the single orbital model, in the degenerated orbital case the spin disorder alone is not 
enough to localize the carriers either, the introduction of the electron-lattice interaction and 
potential fluctuation, within reasonable regions of parameters, can induce the MIT. We also 
discuss the role of orbital disorder in the MIT. It is shown that in the presence of strong 
electron-lattice interaction, electrons become more difficult to be localized in the orbital 
disordered state than in the orbital ordered state due to the split of orbital degeneracy and 
the anistropic electronic hopping in the latter. 

II. MODEL AND METHOD OF CALCULATION 

We consider an effective DE Hamiltonian by taking into account the strong on-site Hund's 
coupling between the itinant Cg electrons and the Mn core spins, together with the coupling 
between electrons and lattice distortion, 

ij i i i 

where c| = (c|+, c|_) with -|- and — representing orbital states ds^'^^r'^ and d^^-y^, respectively. 
The first term of Eq. (1) is the effective DE Hamiltonian in which fij — cos(6'i/2) cos{9j/2) + 
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sin(^j/2) sm{9j /2)e~'^^'^'^~'^^^ with (^j, (pi) being the polar angles characterizing the orientation 
of the localized spin on site i, = are 2x2 hopping matrix in orbital space, with 
a{= x,y,z) being the directions of the bond (ij) between neighboring Mn sites, and — 
t/Aj'^l = &L = 3i/4, = t,i%_ = = 0, and i% = i%_ [13]. The second 

term represents the potential fluctuations experienced by Cg electrons due to the randomness 
of the B?^ and A^^ ion cores, for simplify we assume that is uniformly distributed within 
the range [-W/2, W/2] [8]. The last two terms are related to the Jahn-Teller distortion, in 
which TxiTz are Pauli matrices and sin cos 0^ are the two Jahn-Teller modes [12]. 
Like in many other works, in Eq.(l) we have used the classical limit for both the core 
spins and the Jahn-Teller distortions. In principle, the quantum lattice and spin fluctations 
could change some of the results. In the case of large core spin S — 3/2 and static lattice 
distortion, as is considered in this work, one expects that the quantum effect will not change 
the quanlitative conclusion of the classical approach. 

In the following we will calculate the electronic localization length of Hamiltonian (1) in 
a long bar. Without losing generality we set the bar along the z direction, and the width 
of each slice of the bar is M. Usually the localization length Am is obtained through the 
transfer matrix method [19], which is based on the recursion relationship of successive slices, 
An+i = T~^[{E — Hn)An — T„_iyl„_i], where A^, is the wave function amplitude in the 
nth slice, E is the energy of electron, and is the 2M^ x 2M^ matrix of the electronic 
Hamiltonian inside the slice, and T„ is a 2M^ x 2M^ diagonal matrix with elements being the 
electronic hopping amplitudes between neighboring sites in the nth and n + 1th slices. The 
size of the corresponding transfer matrix is AM'^ xAM'^ . For the present specific Hamiltonian, 
this method can not be directly applied as does not exist due to the lack of hopping 
between neighboring dx2_y2 orbitals along the z direction. The propagation of dx2-y2 in the 
z direction is through its mixing with the o?3z2_^2 in the transverse directions. As a result, 
it is adequate to use the recursion equation of the latter to calculate the localization. By 
eliminating the amplitude of d^^-y^ orbitals, the equation reduces to 
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where is the amphtude of dsz^-r-^ orbital, is the reduced x diagonal matrix 
whose uth diagonal elements is the effective hopping fij between site u (u = 1,2...,M^ 
denote the sites within the slice) in the nth slice and its neighboring in the n + 1th slice. 
H2'^' (7, 7' = +, — ) is the x submatrix of in subspace spanned by basis {cj,^} and 
{cv^/}, Un = H^~{E — H~~)~^H~~^ is due to the mixing of the dsz'^-r^ and d^^-y^ orbitals 
inside the slice. 

The 2M^ X 2M^ transfer matrix of each slice can be obtained from Eq.(2). The localiza- 
tion length corresponds to the smallest positive Lyapunov exponent of the product of the 
transfer matrices, which comes from, if we treat En as a matrix, the behavior of its smallest 
eigenvalue in the limit of n — > 00. In this method, since in a direct computation of the 
product matrix the information associated with the smallest positive Lyapunov exponent 
will be lost when the ratio of the contribution from the smallest positive Lyapunov exponent 
to that of the largest Lyapunov exponent becomes comparable with the machine accuracy, 
one has to orthonormalize the product matrix regularly [19]. Instead of this method, here we 
develop an inverse matrix method in which there is no need to do the orthonormalization. 
Since the smallest eigenvalue of S„ corresponds to the largest eigenvalue of B~^, we seek 
the recursion formula of B~^. Prom Eq.(2), it is given by 

B'l, = 5- V+i (3) 

where the matrix Pn+i can be directly obtained from recursion Pn+i — [E — H;^'^ — Un — 
Tn-iPn\~^T!^. The largest eigenvalue of B^^ is related to the localization length and it in- 
creases exponently with n, so B~^ needs to be regularly (but not necessarily at each n) 
renormalized by divided by hn — \JHij \ {B~^)ij\^. For a bar of length L, the localization 
length is given by A = L/\ YL'n^og{hn)\-, in which the sum is over all n where B~^ are renor- 
malized. The result of A is not sensitive to the initial conditions so that pi and Bi can be 
simply set as the unit matrix. 



6 



The present inverse matrix method has been compared with the transfer matrix method. 
For the usual locahzation problems such as those in Ref.[19], both methods give the same 
results within accuracy. For the present Hamiltonian Eq.(l), we find that the former is more 
efficient due to the smaller matrix size and no need to orthonormalize in the calculation (the 
disadvantage of having to inverse the matrix becomes less serious since in the present problem 
the latter also needs to perform matrix inversion in [/„) . Moreover, it is found that the result 
of the latter lacks stability, for example, the results can be totally different by using single 
or double precision calculations, this deficiency may come from for very small effective 
hopping. On the other hand, there is no such problem in the former method. The difference 
in stability between these two methods is the main reason that we choose the inverse matrix 
method in our calculation. 

III. RESULTS AND DISCUSSION 
A. Pure DE effect 

Now let us first see the case of a pure DE model {g — W — 0) with orbital degeneracy. 
In this case, for a perfect FM phase all electronic states are extended, with the energy band 
from -3t to 3t. For hole-doped manganites Ri_a;Aj;Mn03 (0 < x < 1), the Fermi energy ep 
is at the bottom of the band near x — 1, increases with the decreasing of x, and reaches the 
center of the band at a: = 0. With the increase of the randomness of spin configuration, the 
effective electronic hopping and energy band width decreases. At the PM phase, if we treat 
cos^j and (fi as independent variables that are uniformly distributed in regions [—1, 1] and 
[— 7r,7r], then the band bottom is at around —2.5t. Under this treatment. Fig. 1 shows the 
calculated rescaled localization length as a function of energy E with different width M of 
the bar. The length of the bar is taken to be L = 10^, which is enough for our calculation. 
All curves are crossed at a fixed point Ec ~ — 2.0i, which corresponds to the Fermi energy 
at doping Xc ~ 0.85. The localization length in the thermodynamical limit M ^ oo can be 
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obtained through the one-parameter scahng theory, in which the rescaled locahzation length 
Xm/M can be fitted in terms of a universal function \m{E)/M = f{M/^{E)) [19], and the 
scaling parameter ^ has two branches (inset of Fig. 1). For E < Ec, Xm{E)/M decreases 
toward zero and ^ is the localization length Aoo(-E') for an infinite system, so electrons are 
localized and have no contribution to the conductivity. The electron state becomes extended 
when E > Ec, where Xm/M increases with M and Aoo is infinite. The present calculation 
indicates that in a pure DE model with orbital degeneracy, at doping x < 0.85 in the PM 
state the DE spin disorder is not enough to locahze the carriers on the Fermi surface and the 
system is metallic. This conclusion is not consistent with the experimental results, where 
the MIT usually occurs in the range 0.17 < x < 0.4. 



B. Joint effect 

Next we take the other two scattering mechanisms into account. At each site the JT 
distortion is described by two quantities Qi and 0j. For the amplitudes of the distortion Qi, 
by minimizing the energy their average can be obtained as {Qi) = (1 — x)g. To simplify our 
study in the following calculations we will replace all Qi with {Qi). For the angle 4>i, it is 
closely related to the orbital polarization on site i, and since the configuration of orbitals at 
different sites are disordered at hole concentration 0.17 < x < 0.4 where the FM metallic 
to PM insulator transition occurs, we put 4>i as random variables with uniform distribution 
between — tt and tt. Fig. 2 shows the calculated phase diagram for the PM and perfect FM 
magnetic background at x = 0.2 in the g — W plane, the energy is taken to be the Fermi 
energy at each g and W. Here we are interested in the region between the two MIT lines 
Lp]y[ and Lpy^, which is an insulator in the PM and a metal in the FM case, respectively, 
so that an MIT will occur in this region when the magnetic configuration is changed from 
FM to PM. From Fig. 2, at g = 1.6, to obtain the MIT, the required strength of the on- 
site disorder is 1^ f=:i 3.5t, much smaller than that in Ref.[8] (where W > 12t), and is a 
reasonable magnitude according to the estimation in Ref.[4]. Fig. 2 is the main result of 



8 



this paper, based on it we propose that the FM metaUic to PM insulator transition is due to 
the joint effect of the DE disorder, the JT distortion and the on-site potential fluctuation. 
It is worth mentioning that in the present calculation we have neglected the fluctuation of 
the lattice distortion amplitude Qi, which in principle depends on the electron occupation 
state on each site and has a thermal fluctuation around (Qi) . Since the on-site energy due 
to the lattice distortion at each site is proportional to Qi, this fluctuation will introduce a 
potential randomness in addition to that from 1^, so if considered it may further reduce the 
magnitude of W required for localization. 



C. Effect of orbital configuration 

Since it appears that the on-site potential of orbitals d^z^^r'^ and d^^-yi (being —gQi cos 0j 
and —gQi sin (pi, respectively) is more random in the orbital disordered case than in the 
ordered case and thus the diagonal disorder is stronger in the former, one might expect that 
electrons are easier to be localized in the orbital disordered state than in the orbital ordered 
state. To check this point we also perform calculations in the orbital ordered configurations 
and compare them with that in the orbital disordered case. In Fig. 3 the rescaled localization 
length in both the orbital disordered and C-type antiferromagnetic (AF) orbital ordered cases 
are shown, in the latter case (pi is taken to be ±2/37r in the two sublattices, corresponding 
to the d^x^_r2/d3y2_j.2 type order. It is found that contrary to this expection, electrons 
in the orbital ordered case electrons are actually easier to be localized. Similar results 
are also obtained for the A-type AF, G-type AF and FM orbital ordered configurations. 
This finding is consistent with the occurence of the orbital ordered insulator and orbital 
disordered metal observed in experiments [20]. In fact, for each 0j, one can perform a local 
unitary transformation Ui to transform the orbital state basis from |-|-), |— ) to — 
cos(0i/2)|-F) -Fsin(0j/2)|-) and \-')i = - sin(0j/2)|-F) -h cos(0i/2)|-). The electron-lattice 
interaction becomes diagonal under this transformation, with the on-site energy being —gQi 
for state and gQi for |— respectively, which are independent of the orbital polar angle 
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(pi. At the same time the hopping term in Eq.(l) is transformed to fijU^iijUj, so orbital 
disorder actually does not increase the randomness of the on-site potential fluctuation, it 
only affects that of the electronic hopping, which is already random in the PM phase due 
to the spin disorder. To our understanding, the localization effect of the electron-lattice 
interaction comes from its sphtting the energy degeneracy. In the case of energy degeneracy, 
both orbital states and are important to electronic transport and electron at a 
site can hop to either the |-|-') or the |— ') orbital state on its neighboring site. The electron- 
lattice interaction introduces an energy gap A between these two states. With the increase 
of A, it becomes more difficult for an electron at |-|-') state to move to the neighboring |— ') 
state due to their energy difference, so that electrons are easier to be locahzed in this case. 
To understand that the locahzation effect is even stronger in the orbital ordered case, one 
may note that at large A, electron hopping occurs mainly between neighboring |-|-') states, 
whose hopping amplitude is proportional to cos[(0j — 0q;)/2] cos[(0j — 0a)/2] (a = x,y,z 
labels the (ij) bond direction, ^ = ^2%/ 3 and = 0). In the orbital disordered case this 
amplitude fluctuated from bond to bond and the macroscopic transport are isotropic. When 
orbitals are ordered, however, since orbital ordered state is always accompanied by a global 
anistropy [21], there is one direction along which the hopping amplitude is much smaller than 
that along the other direction(s), for example, in the FM orbital ordered case (pi = 0, one 
can show that whatever the magnitude of ^ is, there is always one direction along which the 
hopping integral is smaller than one quarter of that along another direction. In the special 
case of = TT (corresponding to orbital dx2_y2), along the z direction the hopping even 
vanishes. It is the weaker propagation along this direction that enhances the localization 
effect of the whole system. 

D. Temperature variation 

Finally let us briefly see the temperature induced MIT. In a given system the electron- 
lattice coupling g and potential fluctuation W are fixed while the spin disorder is temperature 
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dependent. For our purpose the temperature dependence of the locahzed spin orientations is 
needed. Here we use the mean-field distribution of the orientations, which satisfies the self- 
consistent equation f{9i,T) = Fexp(3(cos6'i) cos9iTc/T). In Fig. 4 we show the localization 
length as a function of the temperature. For fixed parameters g — 1.7 and W — At, it 
is found that the MIT temperature Tmi ~ 0.9Tc is close to, but not exactly at the Curie 
temperature Tc. Note that here the judgement of M or I phase is whether electrons are 
localized, different from that in some other references where it is based on the sign of dp/dT 
so that the MIT temperature is always at the resistivity peak. Above Tmi and below the 
localization length decreases rapidly with the increase of the temperature due to the sharp 
drop of the magnetization, as seen from Fig. 3, so that the resistivity can still increase in 
this temperature region where transport is dominated by variable-range hopping of electrons 
[5]. 



IV. SUMMARY 

In summary, by using numerical scaling calculations for the localization effect in the 
orbital degenerated manganese system, we have shown that within reasonable range of pa- 
rameters the MIT in manganites can be driven by the joint effect of the DE spin disorder, 
the on-site electronic potential fiuctuation, and the JT distortion. We argue that all these 
three interactions are relevant to the FM metallic to PM insulate transition in colossal 
magnetoresistive manganites. 
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FIGURES 

FIG. 1. Rescaled localization length as a function of the energy E. The inset is the scaling 

parameter ^ as a function of E. The length of the bar is taken to be L = 10°. 

FIG. 2. Phase diagram in the PM and FM magnetic background at a; = 0.2. Lp-^ and Lpy^ 
are MIT lines in the PM and FM states, respectively. Regions below Lpy^ and above L-py^ are 
metal and insulator in both the PM and FM cases, while the region in between is insulator in the 
PM state and metal in the FM state. 

FIG. 3. Comparison of the localization lengths in the orbital disordered (open symbols) and 
C-type AF ordered (solid symbols) configurations in the PM state at 5 = 1.4 and W = At. 

FIG. 4. Temperature dependence of the localization lengths at g = 1.7 and W = At. 
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